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Abstract: A study of Rb and R c is presented in the context of 
a Two Higgs Doublet Model (2HDM) with flavor changing scalar 
currents (FCSC). Implications of the model for the p-parameter 
and for b —>■ s'j are also considered. The experimental data on 
Rb places stringent constraints on the model parameters. The 
configuration of the model needed to account for R b is found to 
be irreconcilable with constraints from b —>■ sj and p. In par- 
ticular, if i?£ xp > Rl M persists then this version of 2HDM will 
be ruled out or require significant modifications. Noting that as- 
pects of the experimental analysis for Rb and R c may be of some 
concern, we also disregard -R^ xp and R e * v and give predictions 
for these using constraints from b — > 57 and p parameter only. 
We emphasize the theoretical and experimental advantages of the 
observable Rb+ C = r(Z — > bb or cc)/T(Z hadrons). We also 
stress the role of Re = T(Z — > hadrons) /r(Z — > £ + £~) in testing 



the Standard Model (SM) despite its dependence on QCD cor- 
rections. Noting that in models with FCNC the amplitude for 
Z — > cc receives a contribution which grows with mf, the im- 
portance and uniqueness of precision Z — > cc measurements for 
constraining flavor changing tc currents is underscored. 



1 Introduction and Summary 



For the past several years precision studies at LEP have been providing 
important confirmation to various aspects of the Standard Model (SM) [p]]. 
A notable exception that has emerged is the decay of Z — > bb. It has long 
been recognized that the Zbb vertex is very sensitive to effects of virtual, 
heavy particles 0. Consequently, a deviation from the prediction of the 
SM could prove to be a significant clue to new physics. It is, therefore, 
clearly important to study Z — > bb in extensions of the SM || and pursue 
the resulting implications. In this paper we study these decays in a class of 
Two-Higgs-Doublet Models (2HDM), called Model III §-0, which present 
a natural mechanism for flavor changing scalar currents (FCSC). 
Our focus is the branching ratio of Z — > bb, i.e. Q 

R b = T[Z - hl) (l) 

F(Z — > hadrons) 

It is worth noting that, since Rb is a ratio between two hadronic rates, most 
of the electroweak (EW) oblique and QCD corrections cancel between nu- 
merator and denominator, making it a uniquely clean and sensitive test of 
the SM. Experiment finds 10: 

Rl xp = .2205 ± .0016 (2) 
whereas the SM prediction is IJ 

R s b M = .2156 (3) 

The difference, of about 3a, is a possible indication of new physics. We note, 
in passing, that the related decay Z — > cc has also been measured albeit with 
appreciably less precision |L| 

R e * p = .1543 ± .0074 (4) 
The SM prediction, on the other hand, is [|l| 

r sm = _ 1724 

Thus i?£ Xp also appears not to be consistent with the SM although the de- 
viation is milder (about 2.3a). It is interesting to note that whereas i?^ xp is 
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larger than Rf* 1 , -R^ xp is less than the SM expectation. Note also that i?^ xp 
quoted above is obtained by holding R c fixed to its SM value [|l|]. 

Our findings are that if we take i?^ xp at its face value then, while Model 
III can accommodate Rl* p , the model parameters get severely constrained. 
In particular, the resulting configuration of the model cannot be reconciled 
with the constraints from the p-parameter and Br(B — > X s ~f). 

Several aspects of the Rb, R c experimental analysis are, though, of con- 
cern. The results given above in eqs. (H) and (f|), include systematic errors 
and emerge from combining the numbers from the four LEP detectors [p]]. 
Since some of the assumptions are common, treatment of the systematics 
can be problematic. Also the errors for Rb and R c are correlated Indeed 
i?£ xp + Rf p is consistent with the SM accentuating the possibility that part 
of the effect may well be resulting from misidentification of the flavors. In 
addition, the observable Ri, 

_ T(Z -> hadrons) 

Rl = v{z - in-) (6) 

which is measured much more precisely than Rb or R c and can be predicted 
in the SM using a s {Mz) deduced from other methods (e.g. lattice and/or 
event shapes in e + e~ annihilation), is found not to be inconsistent with the 
SM, at present. 

In light of these reservations we also fix the parameter space by using only 
the p-parameter and Br{B — > X s ^) and predict Rb, R c and Rb+ C in Model 
III. In particular, in this model, with constraints from the p-parameter and 
Br{B — > X s j), we find that Rb cannot exceed Rl M . Thus, if the current 
trend in the experimental numbers (i.e. i?^ xp > Rl M ) persists, this class of 
2HDM will be either entirely ruled out or require a significant alteration. 

In passing we also emphasize the advantages of the observable Rb+ C 

= r(Z -> 66 or cc) _ 
c F(Z — > hadrons) 
and give the predictions from Model III for Rb+ C - 

Finally, we stress the importance of precision determinations of Z — >• cc 
(i.e. R c ). In type III models its amplitude receives a contribution which grows 
with nij. A precise determination of Z — > cc, thus, constitutes a uniquely 
clean method for constraining the flavor-changing tc vertex that is of crucial 
theoretical concern. 
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2 Two Higgs Doublet Model with Flavor 
Changing Currents 

A mild extension of the SM with one additional scalar SU(2) doublet opens 
up the possibility of FCSC. For this reason, the 2HDM scalar potential is 



usually constrained by an ad hoc discrete symmetry ||11|| , whose only role is 
to protect the model from tree-level FCSC. As a result one gets the so called 
Model I and Model II, when up-type and down-type quarks are coupled 
to the same or to two different doublets respectively |jl2| . In particular, it 



is to be stressed that from a purely phenomenological point of view, low 
energy experiments involving K°-K°, B°-B° mixing, Kl — > fifl etc. place 
very stringent constraints only on the existence of those tree level flavor 
changing transitions which directly involve the first family. Indeed, in view 
of the extraordinary mass scale of the top quark, it has been emphasized 
by many that anomalously large flavor- changing (FC) couplings involving 
the third family may exist ||-[10|,fl3]. Thus, following Cheng and Sher 



perhaps a natural way to limit the strength of the FCSC involving the 
first family is to assume that they are proportional to the masses of the 
participating quarks. In this way, the FC couplings are automatically put in 
a hierarchical order and the third family may well play an enhanced role. 

For this type of 2HDM, the Yukawa Lagrangian for the quark fields can 
be taken to have the form [R PI 



A'" 1 >l'jQ,,0:C, H ■ ////(^./.o,/;,/, • C'jQu.o.C, ,; ■ ^Q, + h.c. 

ud (8) 

where fa, for i = 1, 2, are the two scalar doublets of a 2HDM, while and 
£,fj D are the non diagonal coupling matrices. For convenience we can choose 
to express fa an d fa in a suitable basis such that only the r^- D couplings 
generate the fermion masses, i.e. such that 

< * >= ( «/1/5 ) ' <h> =° (9) 
The two doublets are in this case of the form 
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V2 
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V2 



V2H+ 
ff 1 + iH 2 



v + ff u J ' \ i x ° J j 
The scalar Lagrangian in the (if , if 1 , if 2 , iJ ± ) basis is such that[0, 



(10) 



1. the doublet 0i corresponds to the scalar doublet of the SM and if to 
the SM Higgs field (same couplings and no interactions with if 1 and 
if 2 ); 

2. all the new scalar fields belong to the 02 doublet; 

3. both if 1 and if 2 do not have couplings to the gauge bosons of the form 
H^ZZ or H^W+W-. 

However, while H is also the charged scalar mass eigenstate, (if , if 1 , H 2 ) 
are not the neutral mass eigenstates. Let us denote by (H°, h°) and ^4° the 
two scalar plus one pseudoscalar neutral mass eigenstates. They are obtained 
from (H°, H\ H 2 ) as follows 



H° 

h° 
A 



(if — v ) cos a + if 1 sin a 
— (ff — v) sin a + if 1 cos a 



(11) 



H 2 



where a is a mixing angle, such that for a = 0, (if , if 1 , if 2 ) coincide with 
the mass eigenstates. We find more convenient to express if , if 1 and if 2 as 
functions of the mass eigenstates, i.e. 



H° = (if cos a - h° sin a) + v 

ff 1 = (h° cos a + ff° sin a) (12) 

H 2 = A 

In this way we may take advantage of the mentioned properties (1), (2) and 
(3), as far as the calculation of the contribution from new physics goes. In 
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particular, only the 0i doublet and the rjfj and rjfj couplings are involved 
in the generation of the fermion masses, while 02 is responsible for the new 
couplings. 

After the rotation that diagonalizes the mass matrix of the quark fields, 
the FC part of the Yukawa Lagrangian looks like 

45 % = ifjQi^U^R + £?Qi, L <f>2D j)R + h.c. (13) 

where Qi,L, U^r and Dj^r denote now the quark mass eigenstates and 
are the rotated couplings, in general not diagonal. If we define V^r to be 
the rotation matrices acting on the up- and down-type quarks, with left or 
right chirality respectively, then the neutral FC couplings will be 

SSral = { V L ,D )~ l • i U ' D • V R D (14) 

On the other hand for the charged FC couplings we will have 



S charged — S neutral CKM 

^charged = KjKM ' ^neutral (1^) 

where V CKM denotes the Cabibbo-Kobayashi-Maskawa matrix. To the extent 
that the definition of the £^- D couplings is arbitrary, we can take the rotated 
couplings as the original ones. Thus, we will denote by £^- D the new rotated 
couplings in eq. (0), such that the charged couplings in fllBD look like 
fT CKM and^ CKM -e D 

We will assume that the couplings are purely phenomenological pa- 
rameters and compare the region of the parameter space that could accom- 
modate Rl* p with the constraints from other physical processes. For con- 
venience, we parametrize the t^- D couplings in such a way as to make the 
comparison with the other 2HDM easier 



£S' = Ay (16) 

This is very similar to what was proposed and used in ref. @, || || |Toll , but 
we want now to allow the factors to vary over a broad range, constrained 
by phenomenology only. In this way we may be able to see if the experiment 



data lead to some new patterns in the coupling behavior [15 
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3 Implications for and R c 

Let us now focus on the calculation of Rb and R c . The main task is to 
compute the corrections from new physics to the SM Zqq vertex, for q = c, b. 
Suppose the reference SM vertex for a Z — > qq process is 



V-sm • 9w — 



ASM i 1 ~ 7S) A SM i 1 + 7S) 

Q,L 7) I" ^q,R 7) 



qZ» 



(17) 



where c w is the cosine of the Weinberg angle and g w is the weak gauge 
coupling. The presence of new interactions will then modify it into 



V, 



qqZ 



■ 9w - 
-1— 97a. 



A 



q,L~ 



75), A (1 
■ h Ag iR 



75) 



qZ^ 



(18) 



where 



A^^A^+A^ (19) 

is the sum of the original SM contribution plus the new one from the £-type 
scalar couplings. In principle, both SM and Model III radiative corrections 
to the Zqq vertex give origin to one additional form factor, proportional 
to a^ LU q v (the a^ v q v ^ form factor is absent because it would violate CP). 
This magnetic moment-type form factor arises at one-loop and should be 
considered as well. We have calculated it and verified that, as is the case 
in the SM, it is very small, at least three orders of magnitude smaller than 
the leading contributions to A^/rv Therefore, we neglect its effect in the 
following discussion. 

In view of the previous discussion and neglecting all finite quark mass 
effects (m q ~ 0) jTTJ, the generic expression for T(Z — > qq), for q = b, c, can 
then be written as 

nZ - qq) = ^~^M Z ((A,, L ) 2 + (A^) 2 ) (20) 

S W C W 

where all kinds of EW+QCD corrections have been reabsorbed in the redefi- 
nition of the QED fine-structure constant a, of c w (s w ) and of the couplings 
^ q ,L(R)- Moreover, the A q! L(R) couplings contain corrections induced by the 
new FC scalar couplings. 
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In order to compute the corrections to R q from new physics, such as due 
to the scalar fields of Model III, we observe that, since R q is the ratio be- 
tween two hadronic widths, most EW oblique and QCD corrections cancel, 
in the massless limit, between the numerator and the denominator. The 
remaining ones are absorbed in the definition of the renormalized couplings 
a and s w (c w ), up to terms of higher order in the electroweak corrections 
[0, [18|, |19| . As a consequence, the A q i ( fi ) couplings will be as in eq. (|l"8|), 
with Aq^jq given by the tree level SM couplings expressed in terms of the 
renormalized couplings a and s w (c w ). This feature makes the study of Rb 
and R c particularly interesting, because the new FC contributions may be 
easily disentangled in the Zqq-vertex corrections. In fact, the presence of 
new scalar-fermion couplings will affect the W and Z renormalized propaga- 
tors too, giving stringent constraints especially from the corrections to the 
p parameter. However, this is not relevant for the specific calculation of Rb 
and will be discussed in later segments of this paper. 

In light of the preceding remarks, we can express Rb and R c in terms of 
Rf M and R™ as follows: 



1 + 5, 

Kq ~ Kq JiTr 

where 



D _ DSM 1 u 1 /Ol \ 

Kq K « [1 + R™5 b + R™6 c ] K } 



ASM A NEW I ASM ANEW 
r _ ? ^qL^qL + ^qR^qR 

(A^') 2 + (A^) 2 1 ] 

for q — b,c. In eq. (|2~lD, terms of 0((A^^) 2 ) have been neglected and the 
numerical analysis confirms the validity of this approximation. 

In particular, we will have to compute A^^ R ^ and A^J^ in our model. 
In Fig. |I]we show a sample of the Feynman diagrams which correspond to the 
corrections to the Zbb vertex, due to both charged and neutral scalars/pseu- 
doscalars. The Zee case is strictly analogous, up to modifications of the 
external and internal quark states. In our calculation, we will assume that 
the FC couplings involving the first generation are negligible and we will 
consider all the other possible contributions from the new £-type vertices, 
containing both flavor-changing and flavor-diagonal terms (see eqs. (|l3l)- 

We examined all the possible scenarios, varying the scalar masses (Mh, 
Mh, Ma and M c ), the mixing angle (a) and the ^-couplings. The striking 
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Figure 1: Typical corrections to the Zbb vertex due to both charged and 
neutral scalar/pseudoscalar, in Model III. 

result emerging from this analysis is that, in spite of the arbitrariness of 
the new FC couplings, there exists only a very tight window in which the 
corrections from this new physics enhance Rb, to make it compatible with 
the experimental indications. We find maximum enhancement for 

• very large h°bb and A°bb couplings, obtained for 



• the phase a = 0; 

• light and approximately equal neutral scalar and pseudoscalar masses: 
M h ~ Ma ~ 50 GeV (i.e. at the edge of the allowed experimental lower 
bound for M h and M A @); 

• much heavier charged scalar masses, i.e. M c ~ 400 GeV or more. 
Lighter charged masses require even more demanding bounds on the 
previous parameters. 




(23) 
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For these values of the parameters we can get: 



0.2185 <R b < 0.2230 (24) 
i.e. quite consistent with the experimental measurements, i?^ xp = .2219 ± 



.0017 21 



We note that the enhanced coupling (|23"D to the b quark means ~ £j£ 
(with X tt ~ 1). Perhaps this signifies the special role of the third family with 
respect to Higgs interactions. For our purpose, of course, these couplings are 
purely phenomenological. 

The previous set of parameters strictly mimic what was already found in 
the context of Model II, i.e. without tree-level FCNC. Indeed our model can 
be compared to that one when the phase a = 0, and the FC couplings are set 



to zero. In this regime, we confirm the results of Ref. Hl8j |i9fl . The pattern 
of cancellation between neutral and charged contributions is still valid in 
Model III as well. The charged contribution to Aj^yU is negative and tends 
to reduce Rb, while the neutral one, for light scalar masses (M^a < 100 
GeV), is positive and tends to enhance Rb. With an assumption like the 



one in eq. QTBD, the neutral scalar and pseudoscalar vertex corrections are 



suppressed due to their small couplings to the 6-quark, unless Ajf, > 1. Thus 
in order to enforce the cancellation, we have to enhance these couplings as 



in eq. (|23"D as well as to demand the charged scalar to be much heavier than 
the neutral scalar and pseudoscalar. 

The crucial difference between the two models is that Model III, unlike 
Model II, does not provide any relation between £ u - and £ D -type couplings. 
In fact, for ~ 60rrib/v as in eq. fl23|), we have that ^ ~ while in 
Model II would be inversely proportional to £^ and we would have at the 
same time a very enhanced £ D -type coupling and a very suppressed ^-type 
one. This is at the origin of the slightly more demanding bounds we have to 
impose on the parameters of Model III with respect to Model II if we want 
Rb > Rl M - This difference will become even more important in the discussion 
of the other constraints, as we will see in a while. 

Moreover, in Model III there are also FC couplings, such as £^ and 
We note that, as far as Rb is concerned, £^ plays a role only in the charged 
contribution to ^M R ) anc ^' s * nce this contribution is negative, we do not 
want to enhance it. On the other hand, ^ affects both the neutral and the 
charged vertex diagram, thus, in principle, it could play some role. However, 



11 



even with any reasonable enhancement, £j£ does not seem to change the result 
significantly. 

The scenario we find turns out to be greatly modified when we incorpo- 
rate two additional constraints: the correction to the p parameter, and the 
implication for Br(B — > X a j). In fact, in the framework of Model III with 
enhanced coupling, the first one turns out to be very sensitive to a heavy 
M c , while the second imposes a severe restriction on the magnitude of the 
£mJ coupling. Let us illustrate them in turn. 



4 p-Parameter Constraints on Model III 

The relation between M\y and Mz is modified by the presence of new physics 
and the deviation from the SM prediction is usually described by introducing 



the parameter p [20, 22], defined as 



p0 ~ pM co^e w (25) 

where the p parameter reabsorbs all the SM corrections to the gauge boson 
self-energies. We recall that the most important SM corrections at the one- 
loop level are induced by the top quark [TP], 

Within the SM with only one scalar SU(2) doublet pQ ree = 1. In the presence 
of new physics we have 

Po = l + Ap- EW (27) 

where A/?q ew can be written in terms of the new contributions to the W and 
Z self-energies as 

A NEW (f}\ ANEW(n\ 

Ap ° — w (28) 



Using the general analytical expressions in ref. [p3|| , and adapting the dis 



cussion to Model III (making use of the Feynman rules given in Appendix 
0), we find that 
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Ap* EW ~ g-^-j (sin 2 a G(M C , M A , M H ) + cos 2 a G(M C , M A , Af h )) (29) 
where all the terms of order (Mjy Z /M 2 ) have been neglected and we define 



M'jMlj, Ml M\Mlj, M\ 



l H.h lvl H,h A lvl H.h lvl H,h 



The determination of mt from FNAL |24J] allows us to distinguish between po 



and p ~ 1 + p top . From the recent global fits of the electroweak data, which 



include the input for m t from ref. |24| and the new results on Rb, p turns 
out to be very close to unity. For Rb=R e b yip as in eq. (0) and m t = (174 ± 16) 



GeV, ref. [|2|] quotes 



po = 1.0004 ±0.0018 ±0.0018 (31) 

This result clearly imposes stringent limits on the parameters of any ex- 
tended model. In particular, if we refer to Section |3] and evaluate Apg EW for 
the set of parameters which was found to give an enhanced value of Rb, we 
find that 

A/?NEW " ^k^ Ml (32) 

where the neglected terms are suppressed as (M| A /M 2 ) or (M^^/M 2 ). We 
observe that, for a = 0, the contributions of the (f>\ and <p 2 doublets are 
completely decoupled and the new physics contributions come from the 02 
doublet only. The <p\ doublet can indeed be identified with the usual SM 
Higgs doublet and its contribution to po is already included in the SM value 
of p. Using eq. (|3~2"D , eqs. (|27| ) and fl3lj ) lead to the following upper bound 
on the charged scalar mass 

M c < 200 GeV (33) 

The upper bound (|33| ) for M c means that, to retain Rb in the range of eq. 
would require even larger coupling ^ than in eq. (|23|), since the latter 
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was obtained with M c > 400 GeV and since, also, we cannot reduce the 
neutral scalar masses below their experimental bounds. 



5 Implications of b —> S7 

Even more dramatically, the requirement of an enhanced £^ coupling clashes 
with the experimental constraint for BriB — > X s , ~f) 



Br(B -> X sl ) = (2.32 ± 0.51 ± 0.29 ± 0.32) x 10" 4 (34) 

where the first error is statistical and the latter two are systematic errors. 

This is a remarkable difference with respect to other 2HDM, in which 
there is still a small compatibility between an enhancement over Rf" 1 and 
the result for BriB —>■ X s , -f) obtained by the experiment [Tj|. We will not 



consider Model I, because it cannot produce an acceptable answer for Rb, 
since the fermion-scalar couplings in this model are either all simultaneously 
enhanced or simultaneously suppressed. Thus a disparity between neutral 
and charged scalar vertex corrections can never be realized in Model I. In- 
stead, let us focus on Model II and Model III. It is interesting to compare 
what "the enhancement of the coupling" means in these two models. We 
then immediately realize that in Model III this implies a new large contri- 
bution from the neutral scalar and pseudoscalar penguin diagrams and an 
enormous enhancement of the charged scalar penguin diagram, due to the 
link between neutral and charged coupling via eq. ([TBI) . 

To calculate the contribution of ho, Aq and H to the Br(B — > X s j), 
we work in the effective Hamiltonian formalism, thereby including also QCD 



corrections at the leading order |26j . Due to the presence of new effective in- 
teractions, we need to modify both the basis of local operators in the effective 
Hamiltonian and the initial conditions for the evolution of the Wilson coeffi- 
cients. This is a well known procedure for calculating the effect of heavy new 
degrees of freedom which do not appear in the evolution of the coefficients at 
low energy, but only in their initial conditions at an initial scale roughly set 
at /i ~ Mw- We refer to the literature for all the necessary technical details 
I2SI, E? 



In particular, when we include the new heavy degrees of freedom (h 0} A 
and H^), there are two main changes that we need to consider. First, there 
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are now two QED magnetic-type operators with opposite chirality, which we 
denote by Q^ 1 ^ and write as |3(| 



Q? ' L) = ^"W^ 1 ± 75)6^ (35) 

We recall that in the SM as well as in Model II the absence of is a 
consequence of assuming m s /mb ~ 0. In Model III, we do not want to make 
any a priori assumption on the ^-couplings, because of their arbitrariness, 
and therefore both Qj and Qj can contribute to the b — > 37 decay. The 
rate T(b — * 57) will be proportional to the sum of the modulus square of 
their coefficients at a scale fi ~ mj, i.e. 

r(6 - S7 ) ~ (|cf \m b )\ 2 + |Cf } K)| 2 ) (36) 

We observe that, due to their opposite chirality, the two operators Qj R ' L ^ do 
not mix under QCD corrections and, in a first approximation, their evolution 
with the scale /i can be taken to be the same as in the SM (for Q[ R ^) and 
equal for both of them. In so doing, we neglect those operators whose effect 
is sub-leading either because of their chiral structure or because of the heavy 
mass of the scalar boson which generates them. 

The second change concerns the initial conditions for the Wilson coef- 
ficients at a scale fi ~ M\y- C7 (m b ) depend in general on many initial 
conditions. However, for the same reasons explained before, the most rele- 
vant new contributions, due both to neutral and charged scalar fields, mainly 
affect Cj R ' L \M\v) . In the following we will discuss the results of our numeri- 
cal evaluation of both neutral and charged contributions and their impact on 
the decay rate for b — > S7. In particular, we will focus on the rate normalized 
to the QCD corrected semileptonic rate, i.e. on the ratio: 



R = T(B -> X sl ) T(b -> sj) 
T(B — > X c eu e ) T(6 — > cev e ) 

= r^, , F(\cf\ mb )\* + \Cf\ mb )\*) (37) 
irf(m c /m b ) \ > 

where f(m c /m b ) is the phase-space factor for the semileptonic decay and F 
takes into account some 0(a s ) corrections to both B — > X c ev e and B — ► 
X s 7 decays (see ref. for further comments). We also neglect possible 
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deviations from the spectator model prediction of T(B — > X s j) and r(5 — > 
X c eu e ). From eq. ( p7|) a convenient theoretical prediction for BriB — > X s ^) 
can be extracted, to be compared with the experimental result. 

As far as the new FC contributions from neutral scalar and pseudoscalar 
go, they are peculiar to Model III, because they contain FC couplings. Were 
it not for the enhancement of they would be completely negligible. When 
ibb — 60mfe/v however, the h and A penguin diagrams give a sizable con- 
tribution, amounting to about 30% correction to the SM amplitude. This 
is still within the range allowed by the experiments, and constitute a first 
non-negligible point of difference with respect to Model II. 

However, the most striking effect emerges when we consider the charged 
scalar penguin. Let us focus separately on C\ R \M W ) and Cj L \m w ) and 
try to make a direct comparison with Model II. We recall that the charged 
couplings for Model II are given by 



ttm[3U L V CKU M D D R + — *— U R MuV CKU D L 

tan p 



+ h.c. (38) 



V2 

where Mjj and Mo are the diagonal mass matrices for the U-type and D- 
type quarks respectively, and tan/3 = t^/fi is the ratio between the vacuum 
expectation values of the two scalar doublets. The analogous couplings for 
Model III are expressed by eqs. (|13|) and ([15]). 

Both in Model II and in Model III, the new contributions to C\ R \M\y) 
happens to be multiplied by two products of Yukawa couplings, which we 
will denote by ch and Using eq. (||), we derive that, in 

Model II these products of Yukawa couplings are given by 

KMC = ^VZV tt ml 1 



y/2 tan 13 



2 



(Ctf b )T } = ^=-V t :V tb m t m b (39) 

On the other hand, in Model III, using eqs. (0) and ( ]T5|) they can be written 

as 



(cU*cU\{IH) _ T/ * T/ ( tU , (U Ks 
l?st QbJch — v ts v tb I Stt T Set y 
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{&tS)i m = v t *v tb + (es + (40) 

In order to compare the two models, let us use the parameterization intro- 
duced in eq. fllEp and let us set all the FC couplings in Model III to zero, 
namely = and ^ = 0. Then, the couplings in eq. fl4Cf) reduce to the 
following form: 



AC 1 

(Ctf b )i H) = KV*&& = -^V: s V tb {\ tt \ bb )m t m b (41) 



From eqs. ( p9]) and (|1), the different behavior of Model II and Model III 
with respect to an enhancement of the £^-like coupling should be clear. The 
following correspondence holds: 

Model II Model III 



(tu*tu\ . i _ \2 (42) 

KSst <itb)ch ■ tan/3 2 —> *tt 



66 

In Model II, the enhancement of corresponds to the choice of large value 
for tan/?, i.e. to a suppression of the (6it*£ib)ch coupling with respect to the 
(^t*^a)ch one 5 which stays the same, i.e. pretty small. In Model III, on the 
other hand, we just require > 60 to enhance R b , but we do not have any 
reason to reduce X tt , since each coupling is independent and arbitrary. As a 
net result the charged scalar penguin diagram is greatly enhanced in Model 
III, even with ^ = and ^ = 0. If we restate these FC couplings to their 
non-zero value, the situation is even worse. 

Let us now consider C\ {Myy). This coefficient is special to Model III 
since it is normally neglected in Model II in the limit m s /m b ~ 0. It turns 
out to be proportional to the other two possible combinations of Yukawa 
couplings, i.e. 



\^st ^tblch — v ts v tb I y*Sbs TSss I \^,tt ^ tc y tb 
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(&tE)i m = v; s v tb + tf) (& + ^d) (43) 

and constitutes a relevant extra contribution to BriB — > X a j), to the extent 
that the FC couplings, namely £^ and are not negligible. 

From a numerical analysis, we obtain that for M c = 200 GeV, Model III 
contribution is about a factor of 40 larger than the SM amplitude. When M c 
increases to about 3-4 TeV the two contributions become comparable. Thus 
BriB —->■ Xgj) restricts 

M c > 5 TeV (44) 

in this version of Model III with the enhanced coupling of eq. fl2"3p that is 
needed to account for R^. 

Since the Model II prediction for BriB — > X s 7) is already barely com- 
patible with experiment, unless M c is quite big, the previous comparison 
clearly shows that any enhancement of the coupling, i.e. of Rb, cannot be 
accommodated by Model III. 

6 Remarks on the Experimental Aspects of 

Rb and R c \ Rb+ C an d Rt* 

The preceding discussion leads us to conclude that Model III cannot simul- 
taneously satisfy the constraints from the p-parameter, Br(B — > X s ^) and 
i?£ xp . Therefore, the model may well be wrong and/or incomplete. We view 
the model as an illustration of the kind of theoretical scenarios that can result 
from a rather minimal extension of the SM, namely due to the introduction 
of an extra Higgs doublet. The main virtue of the model is that it gives a rea- 
sonably well defined theoretical framework in which experimental constraints 
on flavor-changing-scalar couplings can be systematically categorized. 

While the model may well be wrong, it is perhaps also of some use to 
question the experimental results i.e. i?^ xp (and R e * p ). As alluded to in 
the Introduction, the experimental analysis for Rb and R c are correlated 
|[. The deviation from the SM given in eq. (§) appears quite significant 
(~ 3a), but this is only after the results from all the four LEP detectors, and 
several different data sets are combined, including their systematic errors. 
One interesting aspect of the Rb results is that all the experiments find that 
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Rl xp > Rl M , although the significance of individual data sets is typically ~ 
(1-2) a. The final errors given in eq. (0) include statistical and systematic 
errors. To the extent that the experiments are truly independent, one is 
tempted to interpret that they are confirming each other at least on this 
overall trend. On the other hand, it is also conceivable that this is a reflection 
of the fact that some of the systematics (shared by the experiments) are 
causing the problem. 

Ironically R e b * p and R exp deviate oppositely from the SM values. In fact, 
using ref. fjj we get 



Rl xp + R exp = (.2219 ± .0017) + (.1543 ± .0074) 

= .376 ±.018 (45) 

which is quite consistent with the SM 

R SM + ^SM = ^ ( 4 Q) 

It is then natural to be concerned that the experimental effect could, in part, 
arise from misidentification of flavors. 
Indeed Rb+ C defined as 

R r(Z — > 66 or cc) 
b+c Y(Z — > hadrons) 

is a very useful observable. It shares the theoretical cleanliness of Rb and 
R c : it is insensitive to QCD corrections. It has significant experimental 
advantages, though, as separation between 6 and c (which is often difficult) 
need not be made. As a specific example, when charm or bottom decay semi- 
leptonically, the hardness of the lepton is often used to distinguish bottom 
from charm. With the use of Rb+ C , one only needs to separate these heavy 
flavors from the really light ones (u, d, s). 

Of course R^ c cannot be obtained by adding the existing numbers for 
i?^ xp and i?^ xp and we will have to await a separate experimental analysis for 
that. Meantime, we note that Re given by 

T(Z -> hadrons) 

Ri ~ T(z^en-) (48) 
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for which experimental numbers are available [|1] does contain information 
on T(Z — > bb or cc). Indeed [l|] 

Rf p = 20.788 ± .032 (49) 

is rather precisely known with an accuracy of ~ .15% which is significantly 
better than R^ (0.7%) or R* xp (4.5%). R e , though, does depend on QCD 
corrections. The calculation of Re is outlined in Appendix [B|. 

It is important to observe that, to calculate the SM prediction (i?f M ) we 
need to use a s {Mz) deduced from other physical methods (i.e. not T(Z — > 
hadrons)). In this way, i?| xp can provide another constraint on any global fit 
of the SM. Two independent determinations of a s (Mz), for example, come 



from the lattice p2l EU| and from the event shapes in e + e annihilation [20 



af^Mz) = .110 ±.006 
af e ~(M z ) = .121 ±.006 (50) 

We will use the average of the two: a s (Mz) — .116 ± .006. Using Table in 
Appendix ||, we then get the SM prediction 

R s e M = 20.748 ± .043 (51) 

The error in eq. ([51]) corresponds to the .006 error (to la) estimates on the 
central value of a s (M z ). Comparing eqs. (|^) and (0), we see that Rf M is 
consistent with the experimental number, i.e. within about la of the error 
on the experiment alone. 

In passing we note that if the true a s (Mz) was taken to be 0.110 then 

R([a s {Mz) = 0.110] = 20.706 (52) 

which would start to deviate from the experimental result in eq. ( f29"D at 
the 2.6a level. But, with the current experimental accuracy, this deviation 
only occurs if one attributes essentially no error to the .110 central value of 
a s (M^)[Q. We do not consider it reliable, at present, to reduce the the- 
oretical errors so sharply. It is clearly important, though, that the efforts 
towards improved evaluations of a s (M z ) be continued, as then the experi- 
mental precision on Ri could be used more effectively to signal new physics. 
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Figure 2: Br(B -> X s -f) in Model III. 
and 2a (dot-dashed) is also given. 

7 Disregarding R e ^ v 



The experimental result at la (dashed) 



Given the previous analysis, we want now to reexamine Model III without 
imposing the constraint coming from i?^ xp . Instead, we will give predictions 
for Rf,, R c and Rb+ C from the model, subjecting it only to the p-parameter 
and Br[B -> X s j). 

If we disregard Rl* v , then there is no need to enhance £^ and we can 
take Xbb = 1 in eq. (|16D . In this case, Model III predicts a Br(B — ► X s ^) 
compatible with experiments at the 2cx-level, for M c > 600 GeV, as we can 
see in Fig. ^. As soon as is not enhanced anymore, the contribution of 
the neutral scalars and pseudoscalar is completely negligible. Therefore, both 
the value of the mixing angle a and of the neutral scalar and pseudoscalar 
masses (Mh, Mh and Ma) are irrelevant. In particular, Fig. |2] is obtained 
for a = 7r/4 and values for (Mh, Mh, Ma) resulting from the fit to Apo, as 
we will discuss in a while. Due to the qualitative character of our analysis, at 
this point it sufficies to seek consistency with the experiment at the 2<r-level. 



21 



0.217 r 

0.216 - 



0.215 - 
0.214 




0.211 

0.21 I 1 1 1 1 1 

300 400 500 600 700 800 900 

M c 

Figure 3: in Model III for a = (solid) and a = ir/4 (dashed). The SM 
prediction Rf M = 0.2156 is also given (dot-dashed) for comparison. 



Indeed, we took as reference the SM calculation [ 31] , which is already affected 



by a large uncertainty, and computed only the leading corrections due to the 
new scalar bosons of Model III, i.e. without considering the complete LO 
effective hamiltonian analysis. From Fig. ^| we also note that, for M c > 600 
GeV, Model III is difficult to distinguish from the SM (again within 2a), 
unless the present SM calculation (Br(B -> X sl ) = (1.9 ± 0.6) x 10~ 4 |3J) 
is improved |34| . 

With the requirement of a large M c coming from Br(B — > X s j), we need to 
consider the discussion of po again and modify it accordingly. The charged 
scalar cannot be the heaviest scalar particle anymore, otherwise ApQ EW would 
be as in eq. (|32| ) and would contradict the present global fit result (see eq. 
SlD). As already noted in ref. [19 for Model II, there are two other possible 



scenarios 



M H , M h <M c < M A and M A < M c < M H , M h 



(53) 
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Figure 4: R c in Model III for a = (solid) and a = vr/4 (dashed). The SM 
prediction R^ M = 0.1724 is also given (dot-dashed) for comparison. 



in which A/?q ew , as given by eq. (p9|), turns out to be negative, and has in 
this way the extra advantage of cancelling the effect of the top quark SM 
contribution (see eq. (|26|)). We note that none of the previous scenarios 
would be compatible with an enhanced value of Rb, because in that case Ma 
and Mh would be required to be equal and light (see Section |3|). 

From a direct numerical evaluation of Ap[J EW , we find that there may exist 



many possible sets of mass parameters for which eq. fl3T|) can be satisfied. 
For instance, let us consider the case in which M#, < M c < Ma- The 
other case in eq. ( |53"D has been studied too and it gives comparable results. 
In order to have a small ApQ EW , it is crucial that M c and Ma are not too far 
apart. One possible optimal set of values for the mass parameters is given 
by the following ratios: M H = 0.4 M c , M h = 0.5 M c and M A = 1.1 M c . In 
this case, the results for R b , R c and Rb+ C are illustrated in Fig. [S] - Fig. |5] 
respectively. The SM predictions are also plotted for comparison. Clearly, in 
Model III, Rb is less than i?^ M and R c is larger than i?^ M . Thus, if the current 
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Figure 5: Rb+ C in Model III (solid) compared to the SM prediction (dot- 
dashed). The dependence on a is irrelevant. 



experimental trend for i?^ xp exceeding R% M persists, Model III will be ruled 
out. 



8 Conclusions 

We analyzed the decays Z — > bb and Z — > cc in 2HDM with FCSC, often 
called Model III. We find that i?^ xp places severe constraints on this model. 
It requires that ~ < 60 GeV, with significantly enhanced coupling of 
the neutral scalar and pseudoscalar to bb. This parameter space of the model 
cannot be reconciled with constraints from the p-parameter and Br(B — > 
X sl ). 

Since aspects of the experimental analysis are of some concern, we also 
examined the model by disregarding i?^ xp and we give the predictions for R b , 
R c and Rb+ C in this case. In particular, we find that, if the current trend of 
i?£ xp > i?^ M persists, then this class of models will be ruled out. 

We emphasized the importance of Rb+ C and Ri in our analysis. 
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In view of the fact that in models with FCSC the rate for Z — > cc receives 
a correction which grows with m|, we stress that precise measurements of 
Z — > cc could provide unique constraints on the crucial tc-vertex. 



Acknowledgments 

We acknowledge useful conversations with Louis Lyons, Vivek Sharma and 
Shlomit Tarem. This research was supported in part by U.S. Department 
of Energy contracts DC-AC05-84ER40150 (CEBAF) and DE-AC- 76CH0016 
(BNL). 



References 

[1] Results presented at the International Europhysics Conference on High 
Energy Physics, Brussels, 1995 and at the 1 7th International Symposium 
on Lepton-Photon Interactions, Beijing, China, 1995. See also: LEP 
Electroweak Working Group Report 95-02. 

[2] J. Bernabeu, A. Pich and A. Santamaria, Nucl. Phys. B363, 326 (1991). 

[3] Here is an illustrative sample of some of the existing studies: A. Denner, 
R.J. Guth, W. Hollik and J.H. Kiihn, Zeit. Phys. C51, 695 (1991); A.K. 
Grant, Phys. Rev. D51, 207 (1995); J.T. Liu and D. Ng, Phys. Lett. 
B342, 262 (1995); E. Ma and D. Ng, Phys. Rev. D53, 255 (1996); J.D. 
Wells, C. Kolda and G.L. Kane, Phys. Lett. B338, 219 (1994); J.D. Wells 
and G.L. Kane, preprint |hep-ph / 9 5 1 03 72| ; G. Bhattacharyya, G. Branco 



and W.-S. Hou, |hep-ph/9512239| ; T. Rizzo, talk given at the Heavy 
Flavor Symposium, Haifa (Israel), Dec. 1995; P. Chiappetta, J. Layssac, 
F.M. Renard and C. Verzegnassi, [hep-ph/9601306| ; G. Altarelli, N. di 
Bartolomeo, F. Feruglio, R. Gatto and M.L. Mangano, [hep-ph/9601324 



[4] T.P Cheng and M. Sher, Phys. Rev. D35, 3484 (1987); D44, 1461 

(1991) ; see also Ref. §. 

[5] A. Antaramian, L.J. Hall, and A. Rasin, Phys. Rev. Lett. 69, 1871 

(1992) . 



25 



[6] L.J. Hall and S. Weinberg, Phys. Rev. D48, R979 (1993). 

[7] W.S. Hou, Phys. Lett. B296, 179 (1992); D. Chang, W.S. Hou and W.Y. 
Keung; Phys. Rev. D48, 217 (1993). 

[8] M. Luke and M.J. Savage, Phys. Lett. B307, 387 (1993); M.J. Savage, 
Phys. Lett. B266, 135 (1991). 

[9] D. Atwood, L. Reina and A. Soni, Phys. Rev. D53, R1199 (1996); 

[10] D. Atwood, L. Reina and A. Soni, Phys. Rev. Lett. 75, 3800 (1995). 

[11] S. Glashow and S. Weinberg, Phys. Rev. D15, 1958 (1977). 

[12] For a review see J. Gunion, H. Haber, G. Kane, and S. Dawson, The 
Higgs Hunter's Guide, (Addison- Wesley, New York, 1990). 

[13] T. Han, R.D. Peccei and X. Zhang, Nucl. Phys. B454, 527 (1995). 

[14] CD. Froggatt, R.G.Moorhouse and I.G. Knowles, Nucl. Phys. B386, 63 
(1992). 



[15] A detailed phenomenological analysis jn| shows that for the first family 
one needs \j ~ 0.1 in eq. (|16| ) in order to satisfy the bounds from K°-K° 
and B°-B° mixing over a large region of the parameter space of Model 
III. Therefore, in this work we will assume that the FC couplings which 
involves the first family are negligibly small. 

[16] D. Atwood, L. Reina and A. Soni, in preparation. 

[17] In our calculation all the quarks except the top quark are taken to be 
massless, keeping m q ^ only in the £^' D couplings as in eq. fll6D and in 
the determination of some corrections to Re, when they become relevant 
to the degree of accuracy that we seek. 

[18] A. Denner, R.J. Guth, W. Hollik and J.H. Kuhn, Zeit. Phys. C51, 695 
(1991). 

[19] A.K. Grant, Phys. Rev. D51, 207 (1995). 

[20] Review of Particles Properties, Phys. Rev. D50, 1173 (1994). 



26 



[21] In obtaining the result in eq. (B|]) we had calculated both Z — > bb and 
Z — > cc in Model III. Therefore the corresponding experimental number 
differs a little bit from that in eq. (Q) which is determined by holding 
R c fixed to its experimental value (see ref. [0). 

[22] P. Langacker, |hep-ph/ 941236 1| , to be published in "Precision Tests of 
the Standard Electroweak Model" , ed. by P. Langacker. 

[23] S. Bertolini, Nucl. Phys. B272, 77 (1986). 

[24] F. Abe et al, [CDF], Phys. Rev. Lett. 74, 2626 (1995); S. Abachi et al, 
[£>0], Phys. Rev. Lett. 74, 2632 (1995). 

[25] R. Ammar et al, (CLEO), Phys. Rev. Lett. 71, 674 (1993); M.S. Alam 
et al. [CLEO], Phys. Rev. Lett. 74, 2885 (1995). 

[26] Note that we are considering the case a = 0, as required by the best fit 
of R b . 

[27] B. Grinstein, R. Springer, and M. Wise, Phys. Lett. B202, 132 (1988); 
Nucl. Phys. B339, 269 (1990); R. Grigjanis, P.J. O'Donnel, M. Suther- 
land and H. Navelet, Phys. Lett. B213, 355 (1988); Phys. Lett. B223, 
239 (1989); Phys. Lett. B237, 252 (1990); G. Cella, G. Curci, G. Ric- 
ciardi and A. Vicere, Phys. Lett. B248, 181 (1990), Phys. Lett. B325, 
227 (1994); M. Misiak, Nuc. Phys B393, 23 (1993); K. Adel and Y.P 
Yao, Phys. Rev. D49, 4945 (1994). 

[28] M. Ciuchini, E. Franco, G. Martinelli, L. Reina and L. Silvestrini, Phys. 
Lett. B316, 127 (1993); Nucl. Phys. B421, 41 (1994). 

A.J. Buras, M. Misiak, M. Miinz and S. Pokorski, Nucl. Phys. B424, 
137 (1994). 



[30] We use the same normalization as in ref. p9[ , in order to check our 
result against the discussion of the problem in Model II, as given in that 
paper. 

[31] M. Ciuchini, E. Franco, G. Martinelli, L. Reina and L. Silvestrini, Phys. 
Lett. B344, 137 (1994). 



27 



[32] A.X. El-Khadra, G. Hockney, A. Kronfeld and P. Mackenzie, Phys. Rev. 
Lett. 69, 729 (1992); C.T.H. Davies, K.Hornbostel, G.P Lepage, A. 
Lidsey, J. Shigemitsu and J. Sloan, Phys. Lett. B345, 42 (1995). 



[33] In this regard see the analysis by M. Shifman ( |hep-ph /95 1 1469| ) , using 



the Z-line shape (rather than Ri) and the value of a s advocated by M. 
Voloshin in Int. J. Mod. Phys. A10, 2865 (1995). It is interesting that 
both these observables (Z-line shape and Rg) lead to similar conclusions 
when a s {Mz) = -HO is used. 



[34] The theoretical prediction of Br(B — > X s ^) from ref. |3l[] includes some 
NLO QCD corrections and the result could change in the future by a 
complete NLO analysis. We could have used in our analysis the fully 
consistent LO result for Br(B — > Xg^), which is a little higher (see ref. 



31, E9|), but this would not modify the qualitative results that we are 



giving. 

[35] A^ EW depend very weakly on m#. The numerical values quoted in Table 
HI are for m H = 200 GeV. 

[36] K.G. Chetyrkin and J.H. Kiihn, Phys. Lett. B248, 359 (1990). 

[37] B.A. Kniehl and J.H. Kiihn, Nucl. Phys. B329, 547 (1990). 

[38] Note that our eq. fl5"7|) differs from eq. (19) of ref. 0. We derive eq. (^) 



from ref. pql , where the results of ref. |37| are improved by absorbing 
large log M|/m^ through the use of the running mass rhb(Mz). Using 
eqs. (2) and (4) of ref. [|36| and substituting the on-shell mass nib with 
the MS mass as in eq. (9) therein, we confirm their result in eq. (14). 
Thus we arrive to our eq. (|57|). 



28 



A Feynman rules for Model III 



In this appendix we summarize the Feynman rules for Model III which are 
used in many of the calculations presented in the paper. 

A.l Fermion- Scalar couplings 

We present the Feynman rules for the couplings of the scalar fields H l (neu- 
tral scalar), H 2 (neutral pseudoscalar) and H + (charged scalar), to up-type 
and down-type quarks, as can be derived from the Yukawa Lagrangian of 
Model III (eqs. (|8|)-([TB|)). Following the discussion of Section these are 
the Feynman rules we need in our calculation of R^. 




-i ((AJ.D . f f/,D*\ . {<-U,D tf [/,D*\ \ 





Although the ^ couplings are left complex in the above, in practice, in 
our calculation we assumed they are real, i.e. t^- D ~ £,fj D * ', as we were not 
concerned with any phase-dependent effects. 
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A. 2 Gauge boson-Scalar couplings 



Here is a list of the Z- and W-boson interactions with Model III scalar fields, 
useful for the computation of ApQ EW . We report them in terms of scalar 
mass eigenstates, H°, h°, A and H + , in order to make contact with the 

We always 



discussion given in Section [| and with the literature |23|, 14 
have to remember the relations (see eqs. (|TTJ) and flI2|)) between the scalar 
mass eigenstates and (H°, H 1 , H 2 , H + ) and use the fact that neither ZH^H 1 
nor ZH°H 2 couplings are present [Ec 



A°,h° 



P2,,' 

it _ ^ ^sina(n 2 -Piy 
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z M T^cosafe 



Pi 



P2,- 



H + ,X + 
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Pi \ ff, x" 
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Pi ^ H\ x" 



30 



H", X 

2 U^2 PI) 



P2 

w5 ^cosafe, -pi)' 



„ H", h° 

P2,'' 



31 



v e, fi,r u,c d,s b 
3.739 2.736 2.200 2.778 -13.848 



Table 1: Values of A^ w , for different flavors, in units of (10 3 ). They have 
been evaluated for m t = 176 GeV and mu = 200 GeV. 

B Calculation of as a function of a s 

In this Appendix we will use the value of a s (Mz) deduced from physics other 
than the width for Z — > hadrons to predict T SM (Z — > hadrons) and i?| M to 
0(a^). Mostly, we follow Bernabeu et al. 0, who give expressions for various 
corrections to T(Z —>■ ff), for both quarks and leptons. 
Let us rewrite the expression for the width of Z — > // as 

T(Z -> ff) = r£(l + A^ w )(l + A^ CD ) (54) 

where Tq is the tree level expression, in which some effects of the EW cor- 
rections have been reabsorbed in the renormalization of the couplings (see 
conventions adopted in 0). A^ w includes only corrections which do not 
depend on a s , i.e. pure EW corrections and QED corrections. They are 
presented in detail in ref. (eqs. (9), (15) and (17), see also references 
therein) and we will not discuss them here. We give their numerical values 
[ |35H in Table |1|. Aq CD represents mostly a s -dependent corrections which can 
be subdivided as 

Aqcd = 5qcd + ^qcd + ^/qcd 

We briefly discuss each of them below. 

The strong corrections to the basic V, A vertex (V = 7 M , A = j^j 5 ) are 
flavor- independent and at 0(af) are given by 

iQCD .^M + 1 . 41 (^M) 2 (56) 

7T \ 7T J 

This is the dominant effect amounting to about 3-4% (see Table 0). 
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SIqcd represents corrections due to kinematic effects of external masses, 
including mass-dependent QCD corrections |36], ^37|. We decide to include 
in the same factor also non QCD mass- dependent corrections to the axial 
vector couplings, in order to make the presentation more compact. Strictly 
speaking, this correction should be included in A^ w . Based on the results 
given in ref. 



36, 371, we deduce 38 



where fi 2 
and 



S f 



3/x 2 



v 2 t + a 2 . 



1 2 

— a* 

2 f 



11 Os 
3 7T 



(57) 



Araj(Mz)/rn 2 y, rrif(Mz) being the running mass at the Z-scale, 



V e 


= -1 + 4x w , 


a e 


= +1 


V u 


= +1 j 




= -1 


Vd 


= 1 + g^W 


ad 


= +1 



Using eq. (2) from ref. ||, we obtain x x 
Numerically, 5 



AtQCD 



-5 x 10~ 3 and 5' 



mQCD 



(58) 



.2314 (where x w = sin 2 ^). 
-0.5 x 1(T 3 (see Table | for 



their a s -dependence). This kind of correction is also relevant, without 0(a s ) 
terms, for the r lepton, in which case it amounts to 5 T ~ —2 x 10 -3 . 

At 0(a 2 ) the large mass splitting between the t and b quarks gives rise to 
a correction, <5/q CD , due to triangular quark loops affecting the axial vector 
current [ 3"7 |: 



°tQCD 



a t CLf 



v 2 + a) 



fOh) 



(59) 



where /(//*) can be written as |}7|> @ 



f(u t ) = log -5 - 3.083 + 0.346^ + 0.211-r (60) 

For m t = 176 GeV we use f{nt) = —4.374. Thus, this correction effects 
+2/3 charge-quarks positively and —1/3 charge-quarks negatively and for 
each flavor it is about 0.4-0.5%, as we can read from Table @. 

Having identified all the corrections to Tf = T(Z — > //), for both quarks 
and leptons, we then consider Re and define 
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a s {M z ) 


$qcd 


°mQcd 


°mQcd 


°iQCD 


°tQCD 


0.105 


34.998 


-5.417 


-0.560 


4.260 


-3.305 


0.110 


36.742 


-5.179 


-0.514 


4.676 


-3.628 


0.115 


38.495 


-4.938 


-0.467 


5.111 


-3.965 


0.120 


40.254 


-4.695 


-0.420 


5.565 


-4.317 


0.125 


42.021 


-4.450 


-0.372 


6.038 


-4.684 



Table 2: Values of different QCD corrections (in units of 10 3 ), for different 
values of a s {Mz)- 



D (T u + T d + T s + T c + T b ) 

Hi = =; (61) 

J- i 

- ST nf ( 1 + A ew) n , A .f \ 

-4- fi T /3) QCD/ 

where Ti = (r e + r^ + r T )/3 and A^ w represents the EW corrections common 
to all the lepton species (see Table [I]). We have denoted by R{ the tree level 
ratios for each quark species. They are given by 



To v 2 e + a 2 e 

< 62 > 

and for x w = .2348 they can be estimated to be R% = 3.4209 and Rf = 
4.4101. 

Finally, Aq CD represents the total QCD corrections for each flavor. They 
are deduced from the previous discussion and their numerical values are 
summarized in Table || together with R^, for different values of a s (Mz)- 

Using the values for A^ w given, for each flavor, in Table |I|, Rg can be 
parametrized as follows (a s = a s (M z )) 



n e,o — 

r>d 

n £,0 — 
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a s {M z ) 


A" 

^QCD 


A c 

^QCD 


A d,s 
^QCD 


^QCD 


Re 


0.105 


39.258 


38.698 


31.693 


26.276 


20.6715 


0.110 


41.418 


40.904 


33.114 


27.935 


20.7060 


0.115 


43.606 


43.139 


34.530 


29.592 


20.7410 


0.120 


45.819 


45.399 


35.937 


31.242 


20.7759 


0.125 


48.059 


47.678 


37.337 


32.887 


20.8108 



Table 3: Values of Ri and its QCD corrections (in units of 10 3 ) as functions 
oia s (M z ). 



R e = Rl (1.000219) (2 + A£ CD (c* s ) + A c QCD (a s )) + 

2R% (1.000796) (1 + A^ D (a s )) + (63) 
R d efi (0.984199) (l + A 6 QCD (a s )) 

from where we deduce the values reported in Table |3|. 
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